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<^ ■ Abstract 

u' 

In this paper we provide a unified approach to a family of integrals of Mellin-Barnes type 

using distribution theory and Fourier transforms. Interesting features arise in many of the 

cases which call for the application of pull-backs of distributions via smooth submersive maps 

defined by Hormander. We derive by this method the integrals of Hecke and Sonine relating 

to various types of Bessel functions which have found applications in analytic and algebraic 

" number theory. 
ON. 
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o: 1 Introduction: 

The theory of distributions and Fourier transforms has been successfully applied in the theory 
of differential operators to obtain precise asymptotic properties of solutions. In this paper we 
look at some applications of distribution theory in the study of special functions by providing 
a unified approach to a certain class of integrals of a type studied by Mellin and Barnes. 

Integrals involving products of gamma functions along vertical lines were studied first by 
Pincherle in 1888 and an extensive theory was developed by Barnes [2] and Mellin jT6]. Cahen 
[5] employed some of these integrals in the study of the Riemann zeta function and other Dirich- 
let series. In the spirit of Mellin's theory some of Ramanujan's formulas have been generalized 
by G. H. Hardy JTU], p. 98 ff. The work of Pincherle provided impetus for the subsequent 
investigations of Mellin [TB] and Barnes [2] on the integral representations of solutions of gen- 
eralized hypergeometric series (see [17], chapter 16 and the comment on p. 225). A detailed 
commentary on Pincherle's work [19] set against a historical backdrop is available in [15]. 

Among the integrals studied by Barnes, the integral formula (13. 6p is well known which served 
as the point of departure for Barnes for his development of the theory of the hypergeometric 
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functions in his seminal paper [2J. A more exotic example fl3.18p - fl3.19p appeared in a later 
paper by Barnes. The integrals of Mellin-Barnes also play an important role in the theory of 
the q-analogues of the hypergeometric functions introduced by Heine in 1847. For a discussion 
of Mellin-Barnes integrals we refer to [26], p. 286 ff and the thorough investigations in the 
recent books [T] and [18]. For the classical evaluation of these and other integrals of this class 
see [1] particularly pp. 89-91 and pp. 151-154. A complete account of q-hypergeometric series 
and their Mellin-Barnes integral representations is available in [9], chapter 4. 

In his work on real quadratic number fields (see [llj, p. 349) Hecke employs the transfor- 
mation formula ( 14. ip of a Mellin-Barnes integral which was generalized by Rademacher with 
a view towards applications to number fields of higher degree (see |20j, p. 58). The integral 
of Hecke also features in his proof of Hamburger's theorem on the Riemann zeta function (see 
[llj, p. 378) where its relation to the Bessel function of the third kind, the Hankel functions, 
is established. Hecke's formula is reminiscent of an integral ( fj4.7p below) considered by Sonine 
in his researches on the Bessel function. The formulas of Sonine and Hecke combine to yield 
an integral representation of the Macdonald function K p (x) that is occasionally employed in 
analytic number theory. 

We provide a transparent and unified approach to these results in the present paper using 
Fourier integrals with non-linear phase functions that yield distributions (densities) given by 
the pull-back of the Dirac delta distribution. We present in sections 3 and 4 proofs of all the 
formulas stated in the introduction (besides a few others) in the spirit of Fourier analysis. It 
turns out that all these formulas fall out of the basic equation proved in section 2: 

/oo p 
dt I exj)(-itp(ui 7 U2,... J u k ))(j)(u 1 ,...,Uk)du 1 ...duk^27r(p*(5o) 7 (j)}. (1.1) 
-oo JR k 

Here p(u±, U2, • • • , Uk) is a phase function, is an arbitrary member of the Schwartz class iS(IR fc ) 
of rapidly decreasing functions (see [T5J, chapter VII), 8q denotes the standard Dirac delta 
distribution and p*(5o) denotes the pull-back of Dirac delta distribution by p (see [12], p. 136 
or [7J, p. 103). It is useful to write (II. ip in the more suggestive notation 

/oo 
exp (-itp(ui,u 2 , ■ ■ -,Uk)) dt = 2ir5(p(u 1 , ...,u k ) = 0). (1.2) 
-00 

A noteworthy special case of (1.2) is the Fourier-G el 'f and formula [TJ, P- 193: 

expit(u — v) dt = 2tt5(u — v). (1-3) 

The paper is organized as follows. Section 2 contains the proof of the basic formula (II. ip 
that we repeatedly use. Section 3 contains the proofs of the formulas of Ramanujan and their 
generalizations due to Barnes, Mellin and Hardy. Section 4 is devoted to the proof of Hecke's 
formula and its relation to various types of Bessel functions. Finally, a few other examples 
susceptible to the same treatment, are gathered up in theorem (4.7). 
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2 A basic formula from distribution theory: 

We shall prove in this section the requisite formula from the theory of distributions that we 
shall frequently employ We recall first the notion of pull-back of distributions via smooth maps 
defined in [12], pp 134-136. For general background on distribution theory we also recommend 
the recent comprehensive work [7]. 

Theorem 2.1 (Pull-back of distributions by smooth submersions): 

(i) Let X and Y be open sets in M. n and R m respectively and p : X — > Y be a smooth 
submersive map. Then there exists a unique continuous linear map p* : T>'{Y) — > V(X) 
such that p*{u) = u o p for all u G C{Y). 

(ii) In the special case when X = M. n and Y = R the pull-back p*(So) is given by 

p*(8 ) = dS/\Vp\ (2.1) 
where dS is the Euclidean surface measure along the zero locus p(xi,x 2 , ■ ■ ■ , x n ) = 0. 
With this notion we state and prove the main formula we use. 

Theorem 2.2: Suppose that p : M™ — > K is a smooth submersive map and <p(x) is a member 
of the Schwartz class <SfM n ), then 

dt <f)(x) exp(itp(x)) dx = 2n (p* (5 ) , (j)} (2.2) 

Proof: Assume first that has compact support along which J^- > 0. Denoting the integral 
in O by J, 



/oo r 
dt / <fi(x) exp(itp(x) — et 2 ) dx 
-oo JR" 

lim / (j)(x)dx / exp(itp(x) — et 2 ) dt 

JR" J-oo 

v^vrlim^ / <p(x) exp(— (p(x)) 2 /4e) dx 



The assumptions on p justify the change of variable^] 

^ : yi = p(xi,x 2 , ■ ■ ■ ,x n ),y 2 = x 2 , ■ ■ ■ ,y n = x n , 



1 The implicit function theorem gives only a neighborhood on which this holds but we may begin by assuming 
that the support of <ft is contained in this neighborhood and then use a partition of unity. 
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and the first of these defines a function x\ = i[)(yi,x 2 , ■ ■ ■ ,x n ). The integral now transforms 
into 

I = yvrlim^ f (j){^{y ll y 2 ,...,y n ),y 2l ... 1 y n )(^) e~ y * /{4t) dy x dy 2 . . . dy n , 

where J^- is evaluated at V2, ■ ■ ■ , Vn)iV2i ■ ■ ■ , 2/n)- Reseating the variable x/i, 

/ = 2V7r hm / r re yi dyidy 2 . . . dy n 

e-+oJ Mn d 1 p{ip{2y/ey 1 ,y 2 , . . .,y n ),V2, ■ ■ ■ ,Vn) 

f (j)(ip(0,y 2 ,...,y n ),y 2 ,...,y n ) , , , 
= 2vr / - - d2/ 2 dj/ 3 . . . dy n , 

Jvn-i (hP{ip{0, V2, ■ ■ ■ , yn), 2/2, • • • , Vn) 



by the dominated convergence theorem and smoothness of the functions involved. Now the 
Euclidean area measure on the locus p(y±, y 2 , . . . , y n ) = is given by 



dS 

3=2 

and invoking the preceding theorem we get 



^ Qljj 2 1/2 I ^/ I 

+ } dy 2 ...dy n = ^dy 2 ...dy n 



f dS 

I = 2n (fr(i)(0,y 2 , . . . ,y n ),y 2 , . . . ,y n )—— = 27r(p*5 ,<i)). 

Jr«-i \yp\ 

The assumption that J^- > along the support of <p is easily removed by employing a partition 
of unity Finally if does not have compact support one can use a sequence (0 m ) of smooth 
functions with compact supports converging to <fi in 5(R n ) and then 



dt \ <f)(x) exp(itp(x)) dx = lim / dt (f) m (x) exp(itp(x)) dx 

= lim 2vr(p*(,5o), <f> m ) = 2n(p*(S ), 0). 



ex- 



Remark: The author is grateful to one of the referees for suggesting the following alternative 
proof of theorem (2.2) using some of the results in [7] and [§]. The notation p* denotes the 
pushforward (see [TJ, p. 93 ff.) under the map p. From problem (10.24) and the Fourier 
inversion formula (14.31) in [7] it follows that 

s -^V 1 ({0}) = 2vrp^ = P*(^lM), 

where 1r denotes the characteristic function of K and J 7 denotes the Fourier transform. Fur- 
thermore we have, for e C7q°(R), 

( p'OFIr), 0) = ( Ji R , p*0) = ( i K , HpM = I r(j>*M)<% 

Jr 
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Using the explicit formula (10.23) in [7] for the pushforward, 

(2^,0) = f [e-** f Mdydxd* 




p-i(W) ngrad || 
e -*p(v) dy dx d£ 
lp-H{x}) ||gradp(y)|| 

Using now the result of exercise (7.36) in [8], the last integral transforms into 
and we get the desired formula (2.2). 

3 Integrals of Barnes, Mellin, Ramanujan and Hardy: 

For details on the gamma function we refer to [7J (pp. 164-168), [17] . [23] and [26] (chapter 12) 
and record here for convenience the basic results that we shall use. 

a\ ^ (a + 1 



Theorem 3.1: (i) Duplication Formula: If Re a > 0, ^vr T(a) = 2°~ 1 r(~)r( 
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(ii) Euler's reflection formula: r(a)T(l — a) = , a Z. 

sin ira 



(iii) For c > fixed, the function t (->■ |T(c + it)\ decays exponentially fast as t — > ±oo 
namely, 

\V(c + it)\ = |t| c "5exp(-7r|t|/2 + 0(l)), t ->■ ±oo (3.1) 

Proof: We shall briefly indicate a proof of (iii) using the Stirling's formula (see [26J, p. 248 
or g] pp. 244-245): 

logr(c-Ht) = (c + it - -) log(c + it) - it + 0(1), |t| -»■ oo. 

Taking the real part of the above equation and exponentiating we get 

\T(c + it)\ = (c 2 + t 2 )§-3 exp(-ttan _1 (t/c) + 0(1)), |t| ->■ oo 

from which (13. ip follows. A different proof is given on pp. 238-239 of [19] (see also p. 234). 
We recall here the definition of the beta function: 

B(p,q)= f u p ~ 1 (l-u) q ~ 1 du, Rep>0, Re q > 0. (3.2) 
Jo 

The basic relation between the beta and gamma function is given by 

T{p + q)B{p,q) =T{p)T{q), Rep>0, Re q > 0. (3.3) 
We begin with a result due to Binet ([1], p. 136) that we shall frequently use. 



Theorem 3.2 (Binet): When Re p > and Re q > 0, 

B (p,q)= / ; r-r^ • (3.4) 

Proof: The result is obtained by substituting u = e x ^ 2 (e x ^ 2 + e" 1 / 2 )" 1 in the integral (13. 2p 
and using the symmetry B(p, q) = B(q,p). 

Ramanujan gave the following formula for the Fourier transform of \T(a + it)\ 2 for positive 
real values of a which we now state and prove using (1.1) (see also j6]). 

Theorem 3.3 (Ramanujan): If a is a complex number in the open right half plane, then 
the following holds. 

f°° 1 

/ T{a + it)T(a - it)e~ l? *dt = x/^r(a)r(a + -)sech 2a £/2, (3.5) 

J -oo 2 

Proof: Using (13. 2p . (13. 3p and (12. 2 p in succession we get 

roo poo poo it(2u-£) J.. 

/ T{a + it)T(a - ity-^dt = 2r(2a) / dt / - _ = AnT(2a)(p*S , (e u + e~ u ) 2a ), 

J —oo J — oo J — oo r T 6 ) 

where p(u) = 2u — £ and so p*(#o) = 5^/2- Thus the integral in ( 13. 5 p equals 

2vrr(2a) 
(4 a cosh 2a £/2)' 

Applying the duplication formula we get the desired result. 

Theorem 3.4 (Barnes' first integral): For a, b, c and d in the right half plane we have 
1 r° „, *. , _ , , , r(a + c)r(a + rf)r(6 + c)r(6 + c/) , , 

— / r a + it r (6 + it r c - it r d - it dt = ^ ; ) ( 1 '-. 3.6 

27r < /_ oc l(a + o + c + a) 

Proof: Assuming a, 6, c and d are real and positive, we split the four gamma factors into 
groups of two and obtain using Binet 's formula, 

cosh(a — c + 2it)u du 



r(a + it)r(c-it) = T(a + c) J 
T(b + it)T(d-it) = T(b + d) J 



oo {e u + e - u ) a+c 
;xD cosh(6 -d + 2it)v dv 



{e v + e~ v ) 



v\b+d 



From this we arrive at the following 

/oo 
T(a + it)T(b + it)T(c - it)T(d - it) dt = 
-oo 

r(0 + c) r(6 + d) / 1 {2cosh J:^ coshv)M £ m it. 
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Here E(t) is the sum of four terms e(t, u, v) + e(t, u, —v) + e(t, —u, v) + e(t, —u, —v) and e(t) is 
give by e(t) = exp((a — c)u + {b — d)v + it(2u + 2v)). Using (11. ip we conclude 



E(t) dt = e {a - c)u+{h ~ d)v 5{2u + 2v = 0) + e {a - c)u - {b ~ d)v 5{2u - 2v = 0) + 
e- {a - c)u+{b - d)v 5(2u - 2v = 0) + e- (a - c)tt - (f> - d),; <5(2w + 2v = 0). 
So we have finally the integral 

J^r(a + it)T(b + it)r(c-it)r(d-it)dt = 2nT(a + c)T(b + d) ^^^^ 

2vrr(a + c)r(6 + d)r(a + d)r(6 + c) 
r(a + 6 + c + d) 

As a noteworthy special case we obtain by virtue of the duplication formula the following result 
and PU, p. 103: 



Corollary 3.5 (Ramanujan): If a and b are positive real numbers, 

-= / \T{a + it)T{b + it)\ dt = x 3.7 

We now turn to an example where the gamma factor appears in the denominator (see [TU], p. 
103). 

Theorem 3.6 (Ramanujan): Suppose Re (b — a) > |, then 

1 f°° T{a + it)T{a-it) F{a)T{a + i)r(6 - a - §) 



v^i-oo T(b + n)r(b-it) " r(6)r(6-i)r(6- a )) ' l ' 

Proof: To first show that the integrand in (13.81) is in L X (]R). Using (13. ip . 



r(a + it)T(a — it) 



< c\t\ 2Reia - b) exp Utan" 1 ^) -Uan-^i/a)) = 0(|t| 2Re(a " fe) ; 



F(b + it)T(b - it) 
from which the assertion follows. We now write 

T(a + it) B(a + it,b — a) T(a — it) B(a — it,b — a 



T(b + it) F(b-a) ' F(b-it) F(b-a) 

and hence the required integral is given by 

4 f f dudv f°° (cosh(2a - b + it)u) (cosh((2a - b - it)v) 



(T(b - a)) 2 J J R 2 (e u + e~ u ) b (e v + e~ v ) b J _ 00 (e u + e- u ) u (e v + e~ v )- u 



dt 
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The integral with respect to t may be written as 

(E(t, u, v) + E(t, -it, v) + E(t, it, -v) + £(t, —u, -v )) dt, 



where J5(t, u, v) = e (^-b)u+(2a-b)v exp _ v _ i og ( e « + e -«) + i og ( e ^ + e -«)}. Note that the 
phase functions are not linear in u and t> . The zero locus of the phase function p(u, v) appearing 
in E(t, u, v) namely, 

p(u, v)=u-v- log(e" + e~") + log(e" + e~ w ), 

is the line u = v along which | Vp(w, i>) | = 2v / 2e~"/ (e n + e~ u ) whereby we get for the pull-back 
the weighted Lebesgue measure p*(#o) = e u (e u + e~ u )du/2. The contribution to the integral 
from Ei(t) is then 

7T /" e ( 4a - 2fe )"e u du 



{T{b - a)) 2 7 R (e« + e-") 2 ^ 1 

Likewise we determine the contributions from E%(t) , Ea(t) , E±(t) and upon adding we get the 
value of Ramanujan's integral (13. 8 j) as 



(r(6 - a)) 2 (e« + e-") 26 " 1 

Using Binet's expression (13.41) and the duplication formula we get (13.81) . 
A more general result may be proved along the same lines: 

Theorem 3.7: If Re (b - a) > 0, Re (d - c) > 0, Re (b + d) > 1 and Re (b + d - a - c) > 1 
then 

1 f°° r(a + it)T(c - it) ,._r(o + c)r(6 + d - a - c - 1) 

— tT/7 7 VrTTr , , T (.""^J 



2tt r(6 + it)T(d - it) T(b + d- l)T(b - a)T(d - c) 

We now state the Fourier transform of T(a + zt) /T(& + it) from which we can get another proof 
of the previous result by appealing to Parseval's formula. 

Theorem 3.8: If Re a > and Re (b - a) > then T(a + it)/T(b + it) e L P (R) where 
{Re (b — a)}" 1 < p < oo. The Fourier transform is the function in L q where q^ 1 > 1 — Re (b — a) 
and is given by 

1 r° T(a + it) - m ,. f e^(l-e^) 6 — \T(b-a))-\ £<0 

^J^WT^ 6 dt -\o. £>o (3 ' 10) 

Remark: Integration by parts reveals that, for £ ^ 0, the integral in (13. lOj) is convergent 
(though not always absolutely) using the following result (see [26], p. 251): 

T'(z) , 1 r 00 wdw 

logz- — -2/ t-it—^TZ^, — tt» Re2:>0. 



r(z) to 2^ Jo ( M 2 + 2 2 )(e 2 ^-l) 

Integral (13.101) is equivalent to formula (13.121) obtained by G. H. Hardy (see [10], p. 98) which 
we discuss next. Integral (13.121) generalizes a formula of Ramanujan. 
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Theorem 3.9 (Hardy): (i) If — p < a < q — p, the following holds: 

j_ r y H«+it )T{a+p+it)T{q _ p _ a _ it)dt =n<^ {31V 



2tt 

ii) If a > and < x < 1 then, 

i rn<>+*)*—*_ ri _ xr w 9 ) (3.i2) 



2n J_ 00 T(a + q + it) 
If x > 1 the value of the integral is zero. 

Proof: (i) We readily transform the integral in ( I3.1ip in the form (1.2): 

POO 

e -u-v u p+a-l vq -p-a-l dudv / exp(^p(«,u)) di. 



y— a rx r x 

~2T 



./0 



The phase p(u, v) = log u — logv — logy vanishes along the line u = vy whereby p*(S ) = vydv 
and the integral is 

POO 

yP / e -v{l+y) v q dv = yPY{q)/{\ + y)1 







as desired. 

(ii) Expressing the integrand in (I3.12p in terms of the beta function and using Binet's formula 
(13.4)) . the integral transforms to 



f 7 — \ , f exp(-it(logx + log(e" + e~ u )) 



2vrr(g) 



^(a—q)u+itu _|_ ^(q—a)u—itu 



dt 



Denoting by p(it) either of the phase functions logo; =F m + log(e M + e u ), we have 

|Vp(u)| = 2|1 -x| 

along the locus p(u) = which is simply the point given by exp(±w) = x/(l — x). With the 
Dirac measure concentrated at this point with normalization (1 — x)/2 we get for the integral 
the value (1 — x) 9-1 /T(q) as asserted. When x > 1 the locus p(u) = is empty which explains 
why the integral vanishes in this case. 

Mellin— Barnes integrals for the Kummer functions: The Kummer functions iF\[a, b, x] 
and U (a, b, x) are normalized linearly independent solutions of the confluent hypergeometric 
differential equation and have the following integral representation (see [22] pp 34-38): 

iF x [a,6,x] = — F ^ I e^t^il-tf-^dt, Re b > Re a > 0, (3.13) 



OC 



U(a,b,x) = p7-T / e~ !Bt * a '~ 1 (l + t) 6 " a_1 ^, Re a > 0, Re x > 0. (3.14) 

1 I (2 ) 
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Theorem 3.10: For real positive values of x we have the following Mellin-Barnes integral 
representation for iFi[a, b, x] and U(a,b,x) 

1 /"» T(b)T(l - it)T(a - 1 + it) 



iFxk^-x] = — / v ') ' \ ^ ... J -x- l+u dt } Re b > Re a. (3.15) 

2* J-oo T(a)T(b-l + tt) 

x i- a r°°r(l-it)T(a-l + it)T(a-b + it) _ it 



U(a,b,x) = — — / ') — - J -x~ lt dt, (3.16) 

K J ^ J-oo r( a )r(i + a - b) y 1 



where Re (a - b + 1) > in fl3~T6|) 



Proof: Expressing the ratio of the gamma factors in the the first integral and using Binet's 
formula we get 

Tib) f°° f°° dv f°° 

' e ' Udu / 1 7, , -„^_i / (E(t, u, v) + E{t, u, -v)) dt, 



27rr(a)r(6 — a)x J J_ OQ (e v + e 



— oo 



where E(t, u, v) = exp(it(logx — logw — log(e" + e ") + v) exp(2a — b — l)v. The zero locus of 
the phase function p(u, v) is the curve x = u(l + e~ 2v ) along which < u < x and 

(/ ' S ' \il + e 2v )du. 



\Vp\ 2 

The sum of the two integrals is then 

x l ~ b T(b) 



e-^\x-uf-^du 

1 {a)L (o — a) J 

from which the result follows. For the case of U(a, b, x), we get 

v ^ b - - T -^r [ e " u " 2iu £ L m - v) + E[t - - v)) dt - 

where E(t,u,v) = exp((a — b — l)v) exp(— it(\ogx — logw — 2v)). The zero locus of the phase 
function p(u,v) = logs — logw — 2v is the curve x = ue 2v along which dS/\Vp\ = \du. From 
this we immediately get the desired result. 

Mellin-Barnes integrals for F(a,b,c, z): Introducing one more gamma factor in f 13 . 1 5 j) 
leads to the following Mellin-Barnes representation for the hypergeometric function (see [2], p. 
142 or [26], p. 286). 

Theorem 3.11: For z not on the negative real axis and Min {Re a, Re b } > 1, the following 
holds: 

F(„ hr 9 \- F ( c ) r T{a-l^it)T{b-\^it)T{l-it)z-^dt 

H a,b,c, z)- T(c-1+U) (3 ' 17) 
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Proof: Taking z = x to be real positive and proceeding as in the case of \F\[a, b, x], the 
integral in (13.17P equals 

I rco rco e - Wvjb -2 dw roo ^ /-oo 

— — — / e u du / — — r-l : — - / (E(t,u,v,w) + E(t,u, —v.w)) dt, 

27rT(b)J J J B(a,c-a)7_ 00 (e- + e-)- i y_ 00 l k ' ; l ' 

where 

E(t, u, v, w) = exp(it(\og(xw) -\ogu + v- \og(e v + e~ v ))e {2a ~ c ~ 1)v . 

The phase function p(u, v, w) is \og(xw) — logw + v — \og(e v + e~ v ) and its zero locus is given 
by xw = u(l + e~ 2v ) which is non-empty only when xw > u. Expressing v in terms of u and w 
we get 

dS xwdudw 
|Vp| 2(xw — u) ' 
and the integral after a scaling transformation assumes the form 

1 

T{b)B{a,c-a) 

— - / 1 5 a - 1 (l-5) C -°- 1 (l+^)- 6 1 



"(1 



"Ms 



-»(l+is) w f'- 



l dw 



Using a familiar result (see [H], p. 245), the last integral equals F(a,b,c, —x). The condition 
that z be real positive may be relaxed by appealing to analytic continuation. 

Theorem 3.12 (Barnes' second integral): Assume that a,b,c,\,fi,u have positive real 
parts and satisfy v = a + b + c + \ + fi. Then the integral 

1 f°° T(a + it)Y{b + it)T(c + it)T(X - it)T(/i - it) dt 

(3.18) 



2tt I> + it) 

equals 

r(A + a)r(A + b)V{X + c)r(/i + q)T(// + bjVjfx + c) 

r(i/ - a)I> - b)T{is - c) [ ' j 

Proof: We derive the result by assuming that a, b, c, A, /i, ^ are real, positive and p > A. These 
may be relaxed later by appealing to analytic continuation. Proceeding as before, 



P( a _|_ j£\ 1 /"°° y e (2a-i/+it)v _|_ g-(2a-!/+it)t> 



r(v + it) T(z/ - a) 7_ oc V (e v + e ~ v ) u+it 

poo poo 

T(b + it)T(c + it) = e^-^wl-^w^ 1 ^ d Wl dw 2 

Jo Jo 

poo poo 

T(\-it)T(fi-it) = e-^-^ul- 1 - 11 ^- 1 - 11 d Ul du 2 

Jo Jo 



dv 
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Thus (I3.18P leads to an integral of the form (1.1) with phase function p given by 

p(v, ui, u 2 , Wi, w 2 ) = \og(w 1 w 2 ) - log(ui« 2 ) ± v — log(e" + e~ v ), iti, u 2 , u>i, w 2 > 0. 
The zero locus p(v, Ui, u 2 , Wi, w 2 ) = is empty if Wi%u 2 < U\U 2 . Along this zero-locus 

dS 2v 



-(1 + e v )du\du 2 dwidw 2 . 



\Vp\ 2 

By virtue of (1.1), the integral (13. 181) equals 

1 r e _ Wl _ W2 _ Ul _ U2 w b-v w c-v u a+\-l u a+»-l, lW2 _ UlU2 y-*~l d Wl dw 2 d Ul du 2 , 

r(^ - a) J R 

where R is the region in the positive orthant given by R = W\W 2 — U\U 2 > 0. After expressing 
it in the form 

poo poo poo rwiW2/ui 

/ du\ I dwi / dw 2 / (j)(ui,u 2 ,Wi,w 2 ) du 2 

Jo Jo Jo Jo 

the integral is simplified in five steps. 

(i) Perform a scaling transformation u 2 = W\W 2 zju\ with z as the new variable (0 < z < 1). 

(ii) Integrate with respect to w 2 . 

(iii) Perform a scaling transformation w\ = su\j z with s as the new variable. 

(iv) Integrate with respect to u\. 

(v) Perform the transformation s = (1 — x)/x with x as the new variable (0 < x < 1). 
We then get the following expression 



T(b + X)T(c + p) f l z a +X -l n_ z y-a-l dz f 1 x *+c-\l- X )b+»-\l- X (l- Z ))-b-*dx. 

T{is-a) J J 

The condition u — a + b + c + X + p, and a familiar result on the hypergeometric function (see 
p. 245) yields 



Et^wwriEO^ r ^ _ + e + v _ l _ dz 

(Y(y-a)) 2 J 
The last integral is again expressible in terms of the hypergeometric functional ([14], p. 277): 

( r( B + A)r( t + A)r(c + A)r( t + rir (c + ^)x + A _ c + v + 

V T(z/ — a)T(A + /i) / 

On appealing to the well known Gauss-Euler formula for -F(a, 6, c, 1) we get (13.191) . 



To apply the formula on p. 277 of [T3] we have assumed fi > A. 
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4 Integrals of Hecke and Sonine: 



A very interesting integral transformation was discovered by Hecke in connection with certain 
problems in algebraic number theory and generalized by Rademacher (see [20], p. 58). The 
formula was later used by Hecke in his proof of Hamburger's theorem on the characterization 
of the zeta function through the functional relation. We discuss in this section this integral and 
a closely related integral representation due to Sonine for the Bessel functions of the first kind. 
We also deduce from Sonine's and Hecke's formulas integral representations of the modified 
Bessel functions. 

Hecke's result gives an Mellin-Barnes integral for the Bessel function of the third kind and 
also ties up with the modified Bessel functions. The most comprehensive work on the theory 
of Bessel functions is |25] but a substantial account is available in the book of Lebedev [H] 
where on pages 117-118 we find a careful specification of the domain of validity for the Hankel 
functions. 



Theorem 4.1 (Hecke): For positive real numbers p, q and arbitrary real a and b, 

27r./_ 00 pa + c + it qb+c+ it ~ J o ex V{ s I s ) s ds ^.1) 

where c is any real number such that c + a > and c + b > 0. 

Proof: Writing out T(a + c + it) and T(b + c + it) as integrals we get 

T(a + c + it)T(b + c + it) dt f°° f°° e^-^u^'W+^dudv f°° . . NX , 

exp(itp(u, v)) dt 



P H Jo Jo P Q 

where the phase function p(u,v) is given by p(u,v) = log(w) — log(pg) and p*(S ) along the 
hyperbola uv = pq is given by the weighted measure pqu~ 1 du. Hence 

r r( ,, + c + «)r(t + c + «)* = 6 _„ r ( _ u _ m„ u 

J_ 00 pa+c+tt q b+c+it Jo V UJ 

The change of variables u = ps now gives the desired result. 

Remark: We note that the special case a = b = and p = q = y/x appears in [21], p. 245 in 
the following form: 

1 f c+lo ° T 2 (s)ds r° e~ 2u ^du 

1 -21 , x > 0. 



27TZ J c _ ioo X s J x yfu< ~ 1 

The substitution u + \/u 2 — 1 = t brings it in the form ( 14. lft . 

Taking p = q we get the Mellin-Barnes integral for the Hankel functions Hp{z) (see [25] . 
p. 180 or [II], p. 108 and p. 118): 
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Theorem 4.2: For the Bessel functions of the third kind we have the following Mellin-Barnes 
integral representation 



me 



^H*H**) = ±E r * +T $£ + a) * «>0 (4.2) 
where r > is arbitrary. 

Proof: Setting p = q = x and b — a = A in (14. X)) and denoting a + c by r we get, 

— / — ; ,o ,o-, — = / exp ( - art - - t A_1 dt 

The change of variables t = e v in the right hand side gives (see [H], p. 118): 

— / — . : Q '— \r. — = 2 exp(-2xcoshu) cosIi(Am) du = tTre tAn/2 H{ >{2ix) 

A further change of variables gives the form in which the result appears in [X Xj , p. 378. 

Remark: The integral formula of Hecke continues to hold when p and q are complex with 
positive real parts. We now derive the Mellin-Barnes integral representation for Bessel functions 
of the first kind (see [3j, Volume - II, p. 21). 

Theorem 4.3: If s is real positive, we have the following 

1 [<°r(\+r + it) , 2r _ m 



J » W = frV-. r(i-r -«) ' *■ (4 - 3) 

where 1 > r > —Re A is arbitrary. 

Proof: In Hecke's formula we set p = q = x + iy with x > and y = — s < 0. We get in the 
limit as x — > 0+ 

1 f°° r(r + it)T(r + A + it) 



2tt 



f°° r(r + U)r(r + A + it) , . , , rr m^, ,n 

/ 1 A+2r + 2^ " ex P("^ + ™r) dt = mH^{2\y\). (4.4) 



Taking the complex conjugate of ( 14. 4 p and performing the change of variables t i— )■ — t we get 
(see jUJ, p. 118): 

1 [°° T(r + it)T(r + \ + it) 
2~ix 

Now subtracting ( 14. 5 p from ( I4.4|) we get 

1 Z" 00 r(r + it)T(r + X + it) 



f°° T(r + it)T(r + X + it) , . , rrf2l/ . .. . , 

/ u A + 2r+2tt -exp(Trt-iTrr) dt = -zvri7f ) (2|y|). (4.5) 

J — oo |y I 



. , |A+2r+2it sin(r7r + i,rt) ^ = 7rJA(2|y|). (4.6) 

z,/l J-oo |y| 

The desired result now follows on appealing to Euler's reflection formula. 

An integral similar to Hecke's with complex values of p and q is the following expression for 
Bessel functions of the first kind due to Sonine (see [25], p. 177): 
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Theorem 4.4 (Sonine): For x > and Re a > we have 

Ja( x ) = / ( c + H)" 1 - 1 exp (c + it- 

where c > is arbitrary. 



.r 2 



4(c + zt) 



dt, (4.7) 



Proof: We prove the result assuming that Re a > 1/2 and then appeal to analytic continu- 
ation. We take advantage of the formula for the Fourier transform of the Gaussian and write 
the term exp(— x 2 /4(c + it)) in the form 

/ 7T \-i/2 r°° 

exp(— x 2 /4(c + it)) = ( ) / exp(—(c + it)y 2 — ixy)dy 

\c + itJ J_ OQ 



and the integral is then 

(x/2) a 



2T(a + ±)tt 3 / 2 



T(a + -)(c + ity a ~2dt / exp Uc + it)(l - y 2 ) - ixyj dy 

-co 2 J-OO ^ ' 



Using the formula 

r>00 

-su„.X— 1 



s - A r(A) = / 

</ 

the integral may be rewritten as: 

(x/2) a 



e~ su u A ~ l du (4i 



c?m / u a 2 exp(c(l — u — y ) — ix?/) d?/ / exp(it(l — u — y )) dt. 



2r(a+|)7r3/ 2 7 

The integral with respect to t is an oscillatory integral with phase p(u, y) = u — 1 + y 2 and we 
get using (II. ip 



-^(p*^),^ 2 exp(c-cw-q/ -ixy))] 



computed along it = 1 — y 2 . Note that since u > only the part of the parabolic locus 
p(u,y) = with y e [—1, 1] is relevant. So we get finally 

(a;/2) a ^ 



v^r(a + §) 



(1 — y ) a 2 cos xyc??/, 

-l 



which is the integral form of Bessel function well suited for obtaining its asymptotic behavior 
(see pj, p. 286 or [14J , p. 114). 

Appealing to analytic continuation, the formula of Sonine holds for all complex values of 
z not on the negative real axis. We now examine the relationship between the formulas of 
Sonine and Hecke, namely from Sonine's formula we derive a representation of the modified 
Bessel function as a Hecke type integral (14. ip . Indeed, Sonine's formula immediately leads to 
the following formula for the modified Bessel functions (see [14J, p. 109 or |20j, p. 47): 



x' p U2x) = — f z~ v ~ x exp (z + —) dz, x > 0, 
2m J L V z / 



(4.9) 
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where p > 0, c > and L is the vertical line traced upwards and passing through c. From this 
we infer the following, where K p denotes the second modified Bessel function (see |14j . p. 108): 

K p (z) = "IzM^lM, HZ; (4.10) 
2 sm pir 

K n (z) = lim K p (z), neZ. (4.11) 

Corollary 4.5: For p arbitrary real we have for all real positive x, 

— sinpn f°° ( x\ du \ r , 
Ip(2x) = / exp ( — xu J — — H — / cospt exp (2x cost) at, 

, v 1 Z" 00 / x\ du 

Kp(2x) = — / exp I — -ux 

2 Jo ^ M 



Proof: Assume first that p is real and — 1 < p < 1. We can then deform the line L in ( 14.9ft 
into a Hankel contour consisting of the pair of parallel rays L± = {z : z = t±ie, — oo < t < — x} 
(e > 0) and a circle of radius x centered at the origin. The Hankel contour is traced so as to 
loop the origin in the counter-clockwise direction. In the limit as e — > we get 

. . sinj97r f°° x, du If 11 

Ip(2x) — / exp(— xu ) — — H / cos ptexp{2x cost ) dt, 

7T J l u u p+L 2n J_ n 

from which we get (by (I4.10p ) the stated integral representation for K p (x). The condition 
— 1 < p < 1 can be relaxed by appealing to analytic continuation since the integrals are entire 
functions of p. The result clearly holds for complex values of x with positive real part. Via 
(14.11) we can immediately get a Mellin-Barnes integral for K p (x) (see [20], p. 55). 

Corollary 4.6: For x and p real positive, we have the Mellin-Barnes representation 

T(z)T{z + p) dz 



L % 2Z 

where L is a vertical line in the right half-plane. 

We now gather up a few other integrals whose proofs, by the method described here, reduce 
to complete triviality. 

Theorem 4.7 (miscellaneous integrals): The following formulas hold: 

1 f°° 

(i) Integral of Cahen: — / T(c + it)y- c - lt dt = e~ y , < c, < y. 



2tt 



-oo 
oo 



(ii) Laplace's integral: — / exp(c + it)(c + it) z dt = (T(z)) 1 , c > 0, Re z > 0. 

2vr ,/_ 00 

/oo 
T(a + it)e~ lt ^dt = 2n exp(—e^ + a£). 
-oo 



16 



Remarks: Laplace's integral appears in his famous work on probability theory (see [13], p. 
137). By deforming the contour of Laplace's integral it is easy to deduce Hankel's integral 
representation for r(z) -1 which is valid for all z G C (see [17], pp. 147-148). The Fourier 
transform (iii) appears in [19] . p. 235. Cahen's integral is valid for complex values of y subject 
to the conditions |arg(y)| < n/2 and y ^ 0. 

We have seen that the method used in this paper handles many integrals which belong to 
the class of Meijer G functions |15j . [3] (volume 1). The precise subclass of Meijer G functions 
to which this method would be applicable is a matter for further investigation. 
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